arXiv:1502.04124v2 [hep-th] 9 Aug 2015 


Bootstrapping the Three-Dimensional Snpersymmetric Ising Model 


Nikolay Bobev^, Sheer El-Showk^’^’^, Dalimil Mazac®, and Miguel F. Paulos^ 

^Instituut voor Theoretische Fysica, KU Leuven, Celestijnenlaan 200D, B-3001 Leuven, Belgium 
^ Theory Division, CERN, Geneva, Switzerland 
^Sorbonne Universites, UPMC Univ Paris 06, UMR 7589 LPTHE, F-75005, Paris, France 
‘^CNRS, UMR 7589, LPTHE, 75005, Paris, France and 
^Perimeter Institute for Theoretieal Physics, Waterloo, Canada 
(Dated: August 11, 2015) 


We implement the conformal bootstrap program for three-dimensional CFTs with M = 2 supersym¬ 
metry and find universal constraints on the spectrum of operator dimensions in these theories. By 
studying the bounds on the dimension of the first scalar appearing in the OPE of a chiral and an 
anti-chiral primary, we find a kink at the expected location of the critical three-dimensional N = 2 
Wess-Zumino model, which can be thought of as a supersymmetric analog of the critical Ising model. 
Focusing on this kink, we determine, to high accuracy, the low-lying spectrum of operator dimen¬ 
sions of the theory, as well as the stress-tensor two-point function. We find that the latter is in an 
excellent agreement with an exact computation. 


Introduction. — Conformal field theories with Af = 2 
supersymmetry in 3d are interesting theoretical models 
with rich dynamics. For example, they enjoy a plethora 
of dualities akin to mirror symmetry in two dimensions 
and Seiberg duality in four dimensions [1-4]. Via the 
AdS/CFT correspondence, they provide a window into 
the nonperturbative structure of M-theory, see for exam¬ 
ple [5]. These theories have also found applications in 
some areas of condensed matter physics, such as topo¬ 
logical phases of matter [6] and optical lattices [7]. 

Of particular interest to us will be the Wess-Zumino 
(WZ) model with N = 2 supersymmetry. This is a the¬ 
ory of a complex scalar (j>, and a complex Dirac fermion 
tf, with the Lagrangian 

Cwz = + + \X\‘^\(j)\'^ 

-I- (A(^V'ae“^'0/3 + c.c.). 

In superspace language, this is the theory of a single 
chiral superfield T = cj) + 9ip + ... with superpotential 
yy = T^. Supersymmetry is not spontaneously broken 
since the Witten index of this theory does not vanish. 
The coupling A is relevant and the theory is believed to 
flow to a 3d AC = 2 superconformal field theory (SOFT) 
in the infrared (IR), which we will denote by cWZ, for 
critical WZ model (see [8] for a review). Note that the 
fermion cannot get a mass term because of the U(l) R- 
symmetry. Supersymmetry then guarantees the scalar 
stays massless as well. From (1), it is clear that this is a 
supersymmetric version of the critical Ising model, whose 
Lagrangian differs only by the absence of fermionic terms. 
This SOFT has been argued to arise as the IR fixed point 
of a certain lattice model [9], and also to describe a quan¬ 
tum critical point on the surface of a topological insulator 
[ 10 ]. 

In this letter, we will use the techniques of the confor¬ 
mal bootstrap to find constraints on the space oi Af = 2 
SCFTs in 3d. These methods have already been applied 
to theories with various amounts of supersymmetry in 


four dimensions [11-16], as well as for AC = 1 [17] and 
Af = 8 [18, 19] theories in three dimensions. The anal¬ 
ysis here is a natural continuation of our work in [20] 
which gives a broad-eyed view of SCFTs with four su¬ 
percharges for general spacetime dimension 2 < d < 4. 
In particular, in [20] we show that certain bounds on 
dimensions of operators show kinks, and we argue that 
one such kink should describe the cWZ model. Here we 
perform a detailed study of this kink for the dimension 
of phenomenological interest, d = 3, following a similar 
analysis for the non-supersymmetric Ising model [21, 22]. 
The outcome is a precise evaluation of the spectrum of 
low-lying operators in the cWZ theory. 

Bootstrap preliminaries .—We are interested in analyz¬ 
ing the consequences of crossing symmetry for a four- 
point function of the form ($$$$), where the operator 
$ is a superconformal chiral primary operator, and $ is 
its conjugate. In other words, $ is the lowest component 
of a short superconformal multiplet and its dimension is 
equal to its R-charge, = <7$. For the concrete exam¬ 
ple of the Af = 2 WZ model in (1), $ will be identified 
with the scalar field </> at the IR fixed point. 

The constraints from crossing symmetry are analysed 
in detail in [20]. Here we only outline the most salient fea¬ 
tures. Decomposing the four-point function by perform¬ 
ing an OPE expansion in the three inequivalent channels 
of fusing $ with $ and $ leads to two crossing equations. 
From the $ x $ OPE one finds that the four-point func¬ 
tion decomposes into a sum of superconformal blocks, 
corresponding to an exchanged superconformal pri¬ 
mary of dimension A and spin s, and its superconformal 
descendants. The blocks take the form 


Ga{u,v) = G%{u,v) - 
s^(A — s — 1) 
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FIG. 1. Bound on the dimensions of the leading unprotected 
operator in the x $ OPE at rimax = 9. There can be no 
unitary SCFTs in the white region. A<s> = 2/3 is indicated 
with a red dashed line. The small shaded rectangle at = 
2/3 indicates the field-of-view in Fig. 2. 

where G^(m, f) is the usual non-supersymmetric confor¬ 
mal block as a function of the two conformal invariant 
cross ratios, u = and v = ^ 31 ^ xhe rep- 

^ 13^24 ^ 13^24 

resentation in terms of ordinary blocks corresponds to 
the four allowed conformal primary operators obtained 
by acting with the supercharges. 

Due to the fact that $ is chiral and using R-charge con¬ 
servation, in the $ x $ OPE, we get a single operator per 
superconformal multiplet, and hence the four-point func¬ 
tion decomposes as a sum of ordinary conformal blocks, 
G^(m, u). The chirality condition also imposes that, in 
this channel, the OPE is of the schematic form (omitting 
conformal descendants) 

$ X $ ~ -p -h ... , (3) 

where Q is one of the Poincare supercharges, IF is an 
antichiral field with dimension 2(1 — A$), and <i>^ is chi¬ 
ral with dimension 2A$. For the range of values of A$ 
considered in this note, supersymmetry imposes a gap be¬ 
tween the dimensions of these operators and the higher 
dimension or spin contributions denoted by ‘... ’ in (3). 
Therefore, when checking crossing symmetry, we should 
allow for operators at these precise dimensions. Since we 
do not know their OPE coefficients, we can only specify 
that the operators may be present, not that they must 
be. See [20] for more details. 

The two crossing equations can be written as an infinite 
number of linear equations with positive coefficients. The 
bootstrap analysis consists of checking under which con¬ 
ditions these equations may have a solution, essentially 
by solving a large linear program. This is done using 
a variation of Dantzig’s simplex method implemented in 
Python [22]. Although the full set of constraints is con¬ 
tinuously infinite, we can truncate the problem to a finite, 
discrete subset obtained by considering finite Taylor ex¬ 
pansions in the cross-ratios u, v around some fixed point. 
We parametrize the truncation by the integer Umax, with 


FIG. 2. Bound on the dimension of the leading unprotected 
operator in the <I> x $ OPE close to =2/3 (upper curves), 
and the corresponding central charge Ct of the solution on 
the boundary (lower curves). The numbers in parenthesis 
next to the values of Umax indicate the number of constraints 
imposed to generate the associated curve. 



Aj 

FIG. 3. A closer view of Fig. 2. The shaded rectangles indi¬ 
cate our estimated error for Aj^j,] and Ct- 

the actual number of constraints (terms in the Taylor ex¬ 
pansion) growing essentially quadratically in this param¬ 
eter (and indicated in parenthesis next to the value of 
n-max in our figures). The reader is referred to [20, 22] 
for more details. 

In the following we consider the problem of determin¬ 
ing an upper bound on the dimension of the leading scalar 
superconformal primary in the $ x I* OPE, which we 
denote by [<i)4>]. Notice that supersymmetry does not 
prevent this operator from having a large anomalous di¬ 
mension. To derive such a bound we impose an ever- 
increasing a gap on the dimension of such operators in 
the OPE until crossing symmetry can no longer be satis¬ 
fied. We only consider a finite number of constraints, but 
since each constraint is a necessary condition for crossing 
symmetry, this yields a rigorous (but potentially sub- 
optimal) upper bound on A[, 5 $]. An important point is 
that if we tune the gap so that we sit precisely at the 
boundary of the allowed region (i.e. maximize the gap), 
it is possible to find a unique solution to the truncated 
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FIG. 4. Bound on the dimension of the subleading supercon- 
formal scalar primary, [‘I>$]', in the <l?x$ OPE. The dimension 
is extracted from the solution that maximizes 


constraint equations [21]. This yields a subset of the spec¬ 
trum of operators and OPE coefficients. As we increase 
the number of constraints this data changes as more op¬ 
erators and OPE coefficients are captured in a convergent 
manner [21]. Below we argue that our bounds are sat¬ 
urated by cWZ, and are able to determine its low-lying 
spectrum with great accuracy. 

Results .— In Fig. 1 we show the results from the boot¬ 
strap analysis which lead to an upper bound on 
There are three distinct features, which we refer to as 
kinks, in the curve bounding the allowed region. The 
second and third kink are discussed in more detail in 
[20]. Here we focus on the first kink, since it bears strong 
resemblance to a similar kink associated with the 3d crit¬ 
ical Ising model [22], and present the numerical results 
pertaining to the small box in Fig. 1. 

A close-up of this region is shown in Fig. 2 : the set of 
curves exhibiting an upward peak are the upper bounds 
on for increasing values of Umax (corresponding to 

imposing successively more bootstrap constraints). Quite 
strikingly, we see that there is a sharp kink for A$ ~ 
2/3. This is precisely the (protected) dimension of the 
chiral field $ in cWZ, since the superpotential W = 
fixes A$ =2/3. Hence it is natural to conjecture that 
our bound is actually saturated by cWZ at this point. 
As we increase the number of constraints imposed, the 
bounds converge rapidly towards this point. This is in 
line with previous bootstrap results that indicate that 
kinks in bounds correspond to actual CFTs. 

In Fig. 2 we also display the central charge [24], ex¬ 
tracted from the OPE coefficients in the unique solu¬ 
tion along the conformal dimension bound curve. Here 
we observe a sharp minimum which occurs very close to 
A$ =2/3. Such minima in the central charge have been 
observed before in the context of the non-supersymmetric 
bootstrap, where they were conjectured to signal the exis¬ 
tence of the ordinary Ising model. Encountering the same 
phenomenon here provides further evidence for our con¬ 


TABLE I. Low-lying spectrum of the critical WZ model. 


A* = 1/2-|-7?/2 

1 (exact) 

— 3 — 1/^ 

1.9098(20) 

A[q-4<[.#] = 3 -I- o; 

3.9098(20) 

A[#s]' 

5.3(1) 

Aj. 

5.25(25) 

Ct 

4.3591(20) 


jecture that the values at the kink correspond to the cWZ 
theory. However, it is important to note that the plot in 
Fig. 2 does not correspond to central charge bounds as 
found in e.g. [22], but rather indicates the value of the 
central charge for the unique solution that maximizes 

A[$$]. 

When examined at high resolution (Fig. 3) we see that 
the horizontal positions of the two kinks do not exactly 
coincide, and they do not occur precisely at A$ = 2/3. 
This is not inconsistent as we have only imposed a finite 
number of constraints (in parenthesis in the plot legends 
next to the value of Umax)- It is clear from Fig. 2 that the 
location of the kinks is changing as a function of Umax 
and seems to be converging towards 2/3. As this note is 
an initial foray into the cWZ theory, we take a conser¬ 
vative approach and estimate the conformal dimensions 
from our most rigorous bounds at A$ = 2/3 and then 
heuristically estimate the error from the rate of conver¬ 
gence. 

Our results for the low-lying spectrum are compiled 
in Table I, where we also present the relation between 
conformal dimensions and critical exponents. To sup¬ 
port our error estimates, we provide, in Fig. 3, 4, and 
6 , plots depicting the convergence rate of these different 
quantities. Note that in the case of A[^|,] our analysis 
yields a rigorous upper bound, since we know A$ = 2/3 
exactly; hence the error bars necessarily lie below our 
bound curve. We believe the analogous result is true 
for the central charge Ct- In Table I, we have also in¬ 
cluded the dimension of the non-superconformal primary 
operator [<5^$$], obtained by acting with the four Q su¬ 
percharges on [$$]. This operator is usually called e' 
in the non-supersymmetric Ising model. Supersymmetry 
implies the relation Ajq4 $$] = A|,|$] -|- 2. Finally, [$<&]' 
denotes the second-lowest nontrivial scalar superconfor- 
mal primary and J' the second-lowest spin-1 supercon- 
formal primary appearing in the OPE of 4> and <i>, the 
lowest being the R-current. Since $, [$4>], [$$]' and 
J' are superconformal primaries. Table I implicitly yields 
the dimensions of their descendant conformal primaries. 

Checks .— Unfortunately, there are not many results in 
the literature for the unprotected spectrum of cWZ. The 
only available result is a 1-loop e-expansion calculation 
[9, 25], which leads to the following low-lying spectrum 
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FIG. 5. Charged scalar spectrum in the vicinity of A# = 2/3. (Left:) the first three spin zero operators in the spectrum 
(colored by order of appearance). The dashed lines correspond to 2A#, d — 2Aj> and 2(d — 1) — 2A# with d — 3 (see [20]). 
(Right:) the OPE coefficients for each operator appearing on the left hand plot (with matching colors). Observe the vanishing 
of the OPE coefficient at A# ~ 2/3. The “noisy” operators in the spectrum plots can be seen to have vanishingly small 
OPE coefficients and thus correspond to small numerical artefacts in the solution. 


[26]: 

A$ = —^= 2 + 0{e^) . (4) 

The expression for A$ is 1-loop exact due to supersym¬ 
metry and after setting e = 1, one finds the expected 
result A$ = 2/3. The fact that Aj$$] = 2 in both d = 2 
and d = 4 (see [20]), together with the expression in (4), 
suggests that Aj^^j never strays too far from 2, and in¬ 
deed this is within 5% of our numerical estimate. 

A strong evidence supporting the identification of the 
kink as the cWZ model comes from the value of the 
stress-tensor two-point function Ct- Dividing by the 
value for the free chiral multiplet Cp®' = 6, the boot¬ 
strap prediction is C^^/Cf"’ ~ 0.7265(3). The same 
quantity can be computed exactly using supersymmetric 
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FIG. 6. The spin-1 spectrum for Umax = 10 in the vicinity of 
A$ =2/3. The decoupling of the subleading spin-1 held is 
reminiscent of what happens in the spin-2 sector of the non- 
supersymmetric Ising model [22]. In the inset we show the 
spin-1 spectrum for various values of Umax which all lie well 
within our estimated error bars. 


localization [27-29], with the result C'^'^ ~ 0.7268, 

showing excellent agreement with bootstrap. The ra¬ 
tio is much smaller than the analogous quan¬ 

tity in the non-supersymmetric critical Ising model [22], 
CtI&t^ ci 0.95. 

An important characteristic of cWZ is the decoupling 
of the chiral operator from the spectrum. This is a 
consequence of the cubic superpotential, which implies 
that vanishes in the chiral ring. Hence, if our con¬ 
jecture is correct, the spectrum at the kink must not 
contain this operator. To this end, in Fig. 5, we examine 
the OPE in the $ x $ channel close to the cWZ kink. 
Notice that the OPE coefficient of <1>^ goes to zero pre¬ 
cisely at the kink. This feature may be thought of as the 
ultimate reason for the existence of the kink, by forcing 
a spectrum rearrangement at this point. The decoupling 
of may also provide an explanation for the peculiar 
see-saw behaviour of the bound in Fig. 1. The next op¬ 
erator in (3) is whose dimension decreases with A$ 
(it is given by Aq 2 ^ = 3 — 2A$). It may be that this 
“drags” A|„j> down, explaining why, to the right of the 
kink, the bound decreases as we increase A$. The de¬ 
coupling of implies also that the first operator in the 
$ X $ OPE is which for A$ = 2/3 has dimension 
5/3 = 1-1- A$. This suggests the identification fk = $, 
so that the OPE takes the form 


$ X 4 ) = -I- ... . ( 5 ) 

This is consistent with cWZ which has a single scalar 
chiral primary operator, and provides further evidence 
for our conjecture. 

In Eig. 6, we examine the spectrum of spin-1 opera¬ 
tors in the $ x $ OPE. Besides the existence of a con¬ 
served current (with A = 2), we observe rapid operator 
rearrangements as the kink is approached. It is inter¬ 
esting to compare this with what happens in the non- 
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supersymmetric Ising model [22, 30]. There, one of the 
prominent features is the decoupling of a spin-2 operator 
with dimension ~ 3.5 when approaching the kink from 
the right. Here we see the supersymmetric analog of this 
transition, with a spin-1 superconformal primary with 
dimension A ~ 2.9 decoupling from the right. In su¬ 
persymmetric theories many spin-2 helds, including the 
stress-tensor, are components of spin-1 superconformal 
multiplets so perhaps it is not so surprising that we ob¬ 
serve a decoupling in the spin-1 sector. It would be in¬ 
teresting to check if this decoupling has a d = 2 analog 
coming from a Virasoro null state (as is the case in the 
non-supersymmetric model). 

Conclusions .— In this note, we have initiated the 
study of the cWZ theory using conformal bootstrap 
methods. Our results hinge on the conjecture that this 
theory sits at the kink in our numerical bounds. Bear¬ 
ing this in mind, we have provided the most accurate 
calculation to date of the critical exponents and OPE co¬ 
efficients in the cWZ model. It will certainly be desirable 
to corroborate our analysis by using other methods, like 
the e-expansion or Monte Carlo estimates. It will also be 
interesting to apply the alternative conformal bootstrap 
algorithm of [31] to 3d A/” = 2 theories. 
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